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Abstract. We have considered scenarios for the evolution of close binaries result- 
ing in the formation of semi-detached systems in which a white dwarf can accumulate 
' Chandrasekhar mass by accretion from a main-sequence or subgiant companion with 

M ~ 2M & . These white dwarfs, probably, explode as type la supernovae or collapse 
with formation of neutron stars. We have carried out a population synthesis study for 
these systems and have estimated the occurrence rate of such events in the Galaxy, 



o 

(N 



p ,■ 

■ depending on the parameter of common envelopes, mass-exchange rate in the binary, 

C/3 ! reaction of the main-sequence component on accretion of helium in the intervening phase 

fNJ ' of evolution. We have found that the model occurrence rate of SNela in semi-detached 

systems is ~ 0.2 x 10 -3 yr _1 , i. e., it does not exceed ~ 10% of the observational estimate 
of the Galactic occurrence rate of SNe la. 
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The problem of progenitors of type la Supernovae (SNe la) is not so lved as yet. Re l- 
atively high fraction of so-called peculiar SNela [possibly, up to ~ 40%. lLi et al. (2001)] 



may suggest that the progenitors of SNela may form a non-homogeneous group. Usu- 
Q H i ally, three possibilities are considered: (i) explosion of a white d warf that has accumu- 

lated Chandrase k har m ass Men by ac cretion in a semi-detach ed dWhela.ri ThenL ll 97.4 



Iben fc Tutukov . 19841 ) or detached ( Tutukov Yungelso Munari fc Renzinil. 
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19921 ) binary system; (ii) explosion of a merger product of a pair of white dwarfs 



with (super) Chandras ekhar total mass dTutukov fc Yungelsonl . Il98ll : IWebbinkl . Il984l : 



Iben fc TutukovL Il984f ): (iii) explosion of a (sub) Chandrasekhar mass white dwarf ini- 



^ ' tiated by the detonation of helium in the laye r of ac c reted matter in a semi-detached 

system with helium donor [edge-lit detonation. iLiviiel (jl990l )]. 

In the present paper we apply population synthesis for close binary stars for the 
estimation of possible occurrence rate of SNela and accretion induced collapses (AIC) 
due to accumulation of Chandrasekhar mass by white dwarfs in semi-detached systems 
with hydrogen- helium donors. We perform a detailed analysis of the dependence of 
results on the parameters of calculations and show that in the Galaxy ^ 10% of all 
SNIa may occur in these systems. 

2 Scenarios of evolution to SNIa 

Scenarios of evolution of close binaries to SNIa are sketched (out of scale) in Figs. H 
and |2 (hereafter, scenarios I and II, respectively). In both scenarios initial mass of the 
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primary component is M\q ~ 5 — 11 M@, mass of the secondary is M20 ~ 0.8 — 3.5 M & . 
In scenario I separation of components corresponds to the RLOF by the primary in 
the hydrogen-shell burning stage (case B of mass exchange). Mass loss by the primary 
results in formation of a common envelope. Evolution inside common envelope leads 
to the shrinkage of the orbit. Due to the mass loss, primary component becomes a 
helium star with a mass ~ (0.8 — 2.5) Mq. Helium stars of this mass rang e expand after 



19851 : iBitzaraki et all I2002I ). This results 



formation of the CO-core (jlben &: Tutukov . 
in the second RLOF and accretion of helium-enriched matter by the companion. Mass 
loss by an expanding star with a CO-core and a helium envelope occurs in the thermal 
time scale of th e star (~ 10 5 yr). Applying results of evolutionary computations by 



Iben fc Tutukov! 1)19851 ) one may estimate the mass of the lost envelope 
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As a result of mass loss, helium star becomes a CO white dwarf. Accretion, most 
probably, leads to the complete mixing o f the secondary (see b elow) . The latter becomes 



somewhat helium-enriched: Y ~ 0.4. Hachisu et al.l l)l999i ) dubbed this scenario "a 



helium-rich supersoft X-ray source channel"]. The secondary fills its Roche lobe in the 
main-sequence stage or soon after developing a helium core. If ensuing mass transfer is 
stable, accreting dwarf may accumulate Chandrasekhar mass. 

Scenario II differs from scenario I by larger initial separation of components. Con- 
sequently, the primary component fills its Roche lobe in the AGB and becomes a white 
dwarf avoiding the stage of helium star. 

In both scenarios in the phase of accre tion onto a white dwarf the system may man- 
ifest itself as a source of supersoft X-rays dKvlafis fc Lambl . ll982l : Ivan den Heuvel et all 
Il992l ). 

In scenario I one of the stages of the evolution is the phase of accretion of helium 
from the envelope of the primary remnant (Fig. stage 5). In our previous papers we 
have assumed that, due to high accretion rate that is defined by Eq. ((21), a common 
envelope that is c ompletely lost by the system is formed in this phase. On the other 



hand, for instance, lHachisu et al.l (|l999l ) assumed that common envelope does not form 
and mass exchange is completely conservative. Actually, the possibility of formation of 
common envelope depends on the reaction of the radius of accreting star on mass gain. 
Thus, this reaction becomes a parameter of the problem. It becomes necessary to study 
the dependence of the model rate of SNIa upon behavior of the accretor radius. 
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3 Accretion of helium by a main-sequence star with initial 



mass of 2M 
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The possibility of accretion of helium, i. e., of the matter with molecular weight 
U exceeding that at the surface of the main-sequence accretor, was noticed long ago 
( Stothers fc Simonl 1969). However, evolution of low- and moderate-mass stars accret- 
ing He was not studied as yet. 

We have considered the evolution of a main-sequence star with initial mass Mq = 
2 Mq that accretes AM = 0.5 Mq of helium-rich matter with a constant rate M a . This 
combination of Mq and AM corresponds to a quite typical close binary system with 
initial mass of components 9 and 2 Mq that may be a progenitor of a SN la [see Eq. I[T|) 
and Fig. El below] . 

Computations were carried out for several values of M a (Table^). The choice for the 
set of M a was defined by the possibility of realization of these rates in the evolutionary 
scenarios resulting in SNIa [see Eq. (j2J) and Fig. Ej. We used for the computations 
the e volutionary code designed for the studies of low- mass stars ( Tutukov Fedoroval . 

Initial chemical composition of accretor was X = 0.70, Y = 0.28, Z = 0.02, the 
chemical composition of accreted matter — Y = 0.98, Z = 0.02. 

We have considered two sets of assumptions on the evolution of accreting star. 

An inverse gradient of molecular wei ght has to le a d to the instability and m ixing 
at the border of helium layer (see, e. g., Ulrichl . 19721 ). iKippenhahn et al, (^)80) have 
shown that for a non-rotating spherically-symmetric star the mixing may be considered 
as formation of helium "bubbles" and their diffusion inside hydrogen-rich substratum. 
Sinking bubbles are gradually destroyed and mixed with surrounding matter. In this 
model, diffusion coefficient D is 
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dr 



(3) 



Here k is opacity, ft is the average molecular weight in the region with the gradient of 
molecular weight and we use the standard notation for the rest of physical variables. 
The time scale r for the doublin g of the thickness W of the transitional layer with the 



gradient of molecular weight is (jKippenhahn et al. , 1980) 



r « W 2 /D. 



(4) 



In the first set of computations we have used expressions (jHJ) and (j3J) for determination 
of the variation of the distribution of He in the star in the course of accretion. 

In the second set of computations we ha ye assumed that the mixin g of accreted mat- 
ter occurs much faster than it was found by Kippenhahn et al. 1 1980h . This assumption 
was based on the estimates of the time scales of physical processes that accompany the 
sinking of turbulent eddies. These time scales, at least in ID approximation, are much 
shorter than the time scale of accretion ~ (10 5 — 10 6 )yr (see Appendix). Therefore, we 
have carried out a set of computations under assumption of an "instantaneous" mixing 
of accreted matter. 
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The parameters of evolutionary tracks are listed in Tabled The tracks of accreting 
stars in the Hertzshprung-Russell diagram and variation of their radii are shown in 
Fig. |21 In the upper panel of Fig. El we show (for track no. 2) the variation of the profile 
of helium distribution inside accreting star with time; this plot also partially includes 
the phase of evolution after cessation of accretion. Computations have shown that under 
assumption of diffusion (tracks 1-3) the radii of stars decrease in the course of accretion 
due to the increase of transparency of the matter with accumulation of helium. 

An expansion of the star that is usually typical for main-sequence stars that accrete 
matter at a high rate occurs under assumption of "instantaneous" mixing (tracks 4-6). 
For accretion rate 2 x 10 -5 MqVT , the maximum radius of the star is almost three 
times larger than the initial one; this ratio is lower for lower accretion rates (Table 
The increase of stellar radii is due to the liberation of internal energy with increase of 
the mean molecular weight of their matter. 

After termination of accretion phase, all models irrespective of mixing mode relax 
to almost similar radii in the thermal time scale (Fig. |IJ lower panel) . 

Thus, evolutionary computations still do not allow to decide definitely what is the 
character of the variation of radii of 1 - 3 M & stars that accrete helium at evolutionary 
defined rates. However, it is clear that the assumption about formation of common 
envelopes due to accretion of helium that was applied in our earlier studies is not the 
only option. 



4 Population synthesis for progenitors of SNIa in semi- 
detached binaries 

In our previous papers we have studied the occurrence rate of accumulation of Mch 
by white dwarfs in the systems of different types: cataclysmic variables that are evolving 
in the thermal time scale of the do nor or in the time scale of angular momentum loss via 
magnetically coupled stellar wind ( Yungelson et al. . 1996); in the systems with subgiant 
donors with d egenerate helium cores that evolve in the thermal or n uclear time scale 



of the donor ((Yungelson et al 



19961 : lYungelson fc Livid ll99Sl . l200Ch . In the former 
systems mass of the donors was limited from above by 1.5 Mq, while in the latter by 
2.5 M & . Our estimates of the occurrence rate of potential SN la in semi-detached systems 



did not exceed ~ 10% of the empirical estimate o f the SNIa rate in th e Galaxy that 



is (4 ±2) x 10~ 3 yr~ 1 dCappellaro fc Turattd l200lh . On the other hand. lHachisu et al 



1 199d and iLi van den Heuvell (|l997l ) have found, that the model rate for the 

systems of the same type may be comparable to the empirical one. We should note, 
however, that the starting point for the estimates made in the papers of Hachisu et al. 
and Li and van den Heuvel were the systems composed by white dwarfs with main- 
sequence or subgiant companions. The authors of the aforementioned papers did not 
estimate the rate of formation of these "progenitor systems " by means of populatio n 
synthesis. This may be one of the causes for the discrepancy ( Yungelson Livid 19981 ). 

In the present paper we study the possibility of accumulation of Mqh by white 
dwarfs evolving according to scenarios I and II, using uniform assumptions on the 
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mass loss by the donor and the possibilities of mass gain by white dwarfs. It is evi- 
dent that most part of the matter potential SNela accrete in the thermal time scale 
of the donor, when the hydrogen burns stationary and the rate of its accumulation is 
M s t — 10~ 7 — Mq yr -1 . In addition to binary systems considered by us before, we 



have included into consideration binaries with main-sequence compone nts that satisfy 
the condition of dynamically stab le mass exchange: q = M2/M1 ,$,2.5 (jTutukov et, al 



19821 : iHiellming fc WebbinMll987h . Ifl.2;$g,$,2.5, mass exchange in these systems pro- 



ceeds in the time scale close to the thermal time scale of the donor until masses of com- 
ponents become equal. Accretion rate M a may be then higher than t he maximum rate 
of sta ti onary hydrogen burn ing M st . Following Kato fc Ibenl (^92); K ato fc Hachisul 
( 19941 ): lHachisu et alJ ((1996J) we have assumed that the stationary burning of the hy- 
drogen generates optically thick stellar wind that takes excess of the matter away from 
the system. Specific angular momentum of the wind matter is equal to the specific 
angular momentum of the dwarf. Mass loss in this mode stabilizes mass exchange in 
the system. Optically thick wind exists only if M a 10 -4 M Q yr _1 . For higher M a 
common envelopes should form. 

To improve the accuracy of the computations of the potential occurrence rate of SN la 
in semi-detached systems, we have introduced some modifications concerning efficiency 
of accumulation of the matter by white dwarfs, mass loss from the systems, and behavior 
of the radii of the accre t ing stars into the populat i on sy nthesis code we have used before 
( Yungelson et al. , 199fil : lYungelson fc Livid . Il998l . l200fll ) . Here we mention only that we 
assume that all stars in the Galaxy are born in bi naries and that we appl y the birthrate 
function of binary stars in the form suggested by llben fc Tutukoy (|l984l ) 
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0.2(i(log a )^^f(qo)dq , 



(5) 
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where ao is the initial separation of components, f(qo) - distribution of initial mass 
ratios of components normalized to 1 (oq = M20/M10). Function Q is normalized to 
the formation of one binary system with Mm > 0.8 per yr. The mass of progenitors 
of white dwarfs is (0.8 - 11.4)M (|Then Tutukovl . 1985). It was assumed that star 
formation rate in the Galactic disk was constant for 13 x 10 9 yr . Initial distributions of 
binaries over q$ and log ao were assumed to be flat. The space of initial parameters of 
the binaries was approximated by rectangular grid with steps 0.0125 in logMio, 0.0125 
in log(ao), and 0.025 in q$. 

The number of systems with mass exchange resulting in the accumulation of Mch 
depends, mainly, on the parameter of the common envelopes a ce . The "equation of 
common envelopes" was taken in the form 



(M a0 + M d0 )AM 



a 



M df M af M a0 M d0 



(6) 



where indices and / indicate initial and final values of the masses of donor and accretor 
and their separations and AM is the mass of the lost envelope. The value of a ce is the 

1 This assumption is justified since we consider rather short-living (~ 10 9 yr) stars with M ^ 1.6 Mq 
(see Fig 01. 
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measure of the ratio of the binding energy of the common envelope and the orbital energy 
of the binary system. At present, neither theory nor observations allow to restrain a ce 
with reasonable accuracy. Our choice of the range of a ce use d in the computations (0.5 to 



witn reasonable accuracy, uur cnoice 01 tne range 01 a ce used m tne comp 
5) is j ustified by t he results of the analyses of observations fib eil 1200(1 F 



Nelemans et al 



2000: ISokeTJ . l2002h that provide evidence for a high efficiency of the deposition of energy 
into common envelopes and for the loss of a considerable part of the donor envelope by 
stellar wind prior to RLOF, that is equivalent to high a ce . Of course, it should be noted 
that the usa ge of the same a™ for al l common envelopes is a forced simplification of the 
problem [see Dewi Sz Taurid (|2000h ; iTauris &: Dewil (|200lh where the binding energy of 
the envelopes is evaluated as the function of the evolutionary state of the star]. 

For every semi-detached system with main-sequence donor that was generated by 
population synthesis code, mass loss rate was assumed to be constant and equal to = 
a ifc^KHi where Mkh is the mass loss rate corresponding to the thermal time scale of the 
donor. It is evident that in real systems depends on the mass ratio of components 
(decreases with decreasing mass ratio of components) and on the evolutionary state of 
the donor at the instant of the RLOF. As well, it varies in the course of evolution. The 
value of ath is, thus, an additional parameter of the model . The value of ath was 
varied from 0.2 to 1. The computations were continued until white dwarf accumulated 
Mch or while mass of the donor became equal to the mass of accretor. In the systems 
with subgiant donors mass loss first occurred in the thermal time scale of the donor and 
then - in the time scale determ ined by the rate of growth of the donor helium core (see, 



Yungelson &: Liyiol . 1998). In these systems computations were limited either by 



e. g. 

accumulation of Mch by the white dwarf or by complete loss of the hydrogen envelope 
by the donor. 

For the estimate of the efficiency of helium accumulation by white dwarfs (or their 
er osion due to Nova e expl osions at M a £ 1.5 x 10~ 8 M Q yr _1 ), the grid of models 
of lPrialnik fc Kovetj (|l995h was applied. The loss of matter via thermal flashes in the 
helium surfa ce layer of white dwarf a ccumulated due to hydrogen burning was estimated 
according to Kato Hachisul (|l999h . 

According to the results of computations of accretion (Sec. EJ) we have considered 
two limiting cases of reaction of the main-sequence star - future donor in the system 
with white dwarf - on accretion of helium-enriched matter: completely conservative 
mass exchange with the radius varying according to the mass-radius relation only, and 
conservative mass exchange that is accompanied by a three-fold increase in accretor 
radius. For the sake of comparison with our previous studies we also considered the 
occurrence rate of SNIa assuming that accretion is accompanied by formation of a 
common envelope. 

Results of computations are given in Table |2j For every combination of otce 

and a t h 

we give separately the occurrence rate of possible SN la accumulating Mch according to 
scenario I under three different assumptions on the reaction of the star upon accretion 
(three upper lines). The fourth line for every set of parameters in the Table gives 
the rate of events in the systems evolving through scenario II. Separately, as AIC, are 
listed accumulations of Mch by white dwarfs with initial mass ^ 1.19 M Q that, under 
our assumptions, consist of the oxygen-neon-magnesium mixture. At accretion rates 



6 



exceeding ~ 1CT 8 M yr _1 these dwarfs, p erhaps, do not expl o de as SNIa but, instead, 
collapse with formation of neutron stars ( Nomoto Kondcl 1991 ). "Initial mass of 
the primary component - initial separation of components" diagram for progenitors of 
the systems in which white dwarfs may accumulate Mch is shown in Fig. |S] [for one 
of the most "fertile" model no. 13 in Table [2] . The same figure shows the relation 
between masses of the donors - main-sequence stars or subgiants and masses of white 
dwarfs in the systems in which white dwarfs are able to accumulate Mq^i at the moment 
when the mass transfer starts. Two groups of objects in the lower panel of Fig. are 
systems with main-sequence donors (M^ > 2.5M Q ) and systems with subgiant donors 
{Md = 1.5 — 2.2 Mq). The cause for the "gap" between two groups of objects is formation 
of deep convective envelopes in the stars with M 2.2 Mq immediately after TAMS; 
this results in dynamical mass loss upon RLOF and formation of common envelopes. 

Interestingly, a significant fraction of "presupernovae" white dwarfs is initially not 
very massive (M a ^ 1M ) and have relatively low- mass companions {Md ^ 2M ). 
The circumstance that in these systems the combination of the mass accretion rate and 
the statio nary burning r ate is favor able for accumulati on of Mch is in agreement with 
results of Laneer et al. ( 2000h and Sarna et alJ (|2002h who have carried out detailed 
computations for the main-sequence star - white dwarf systems. 

The dependence of results on the parameters of the population synthesis may be 
summarized as follows. 



The number of the systems with white dwarfs successfully accumulating Men 
depends, mainly, on a ce . For a ce > 1 the number of systems evolving according to 
scenario I diminishes and for a ce = 2 they practically disappear, because efficient 
deposition of energy into common envelopes limits the separation of components 
from below. As a rule, scenario I is dominating: see Fig. [5] in which the border 
between the systems evolving according to scenarios I and II is shown; this line 
only weakly depends on M20. 

Reduction of ath from 1 to 0.5 results in increase of the number of potential 
SNIa. For high ath a significant fraction of the matter is accreted with the rate 
that is sufficiently higher than the rate of the stationary hydrogen burning and 
is, therefore, lost from the system via optically thick wind. For this reason, the 
efficiency of mass accumulation is low. For low ath = 0.2, efficiency of helium 
accumulation is also low, due to the recurrent explosions of Novae; this results in 
reduction of the number of systems with dwarfs that are able to accumulate Mch- 

In scenario I, for a ce > 1 and ath > 0.5 the assumption about formation of a 
common envelope in the course of accretion of the He-rich matter increases the 
number of potential SNIa, since this increases the number of systems in which 
RLOF becomes possible. For a ce < 1, on the contrary, the mergers of components 
reduce the number of potential SN la. 

The effect of a three-fold increase of the radius of the accreting star as compared 
to the case of mass exchange without expansion of accretor, is comparable to 
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the case with formation of the common envelopes. The reason for the similarity 
is in the narrow interval of the major semi-axes of orbits in progenitor systems: 
A log a « 0.3- 1 (Fig. El). 

Comparison with results of lYungelson fc Livid ( 19981 ) shows that taking into ac- 
count systems with initial mass of the secondary higher than 2.5 Mq results in 
increase of the SNIa rate by less than ~ 30%. The increase of mass of the donors 
more strongly influences the rate of AIC's that almost doubles [c.f. Table 1 in 
Yungelson fc Livid il99Sh ]. 



Semi-detached systems with accreting white dwarfs that burn hydrogen station- 
ary are usually identi fi ed w i th one of the subpopul a tions of supersof t X-ra y sources 
( van den Heuvel et al. . 19921 : Rappaport et al. . 19941 : Yungelson et al. . 1996h . In Ta- 
ble [21 for every set of computations we give the number of dwarfs -/Vgss that may be 
observed as sources of supersoft X-rays (for scenarios I and II together). For the life- 
time of stationary sources we have assumed the time of accretion with M > M s t] for 
nonstationary sources the lifetime was assumed equal to the period of time when, after 
Nova explosion, the lu minosity of the dwarf exceeds 4 x 10 37 ergs _1 [see. Eq. (12) in 



Yungelson et alJ (|1996h ] 



In scenario I that is actual for most systems, the matter of the donors has to be 
enriched in He (Y ~ 0.4). Note, helium-enriched stationary sources of supersoft X-ray 
emission are not known as yet (P. Kahabka, private communication). The only observed 
source enriched in He is recurrent Nova U Sc o. The presence of a white dwarf with the 
mass close to Mch is suspected in this system ( Hachisu et allEbodlTho roughg ood et al 



2001). The estimates of mass exchange rate in this system show that explosions at the 
surface of the dw arf in this system do not result in erosion and the dwarf may really 



accumulate Mch (jSarna et al 1 120021 ; 



In the Galaxy, only two of the sources of supersoft X-rays that are identified optically 
have orbital period s of ~ 16 hr and ~ 4 day and may belong to the subtype considered in 
the present paper (|Kahabkal .E002 N ). The number of sources that may be identified very 
strongly depends on t he ill known interstellar absorpt i on an d the distribution of sources 
over the temperature. Di Stefano Sz Rappaportl (|l994l . Il995l ^l : IPi Stefano fe Kond {2003) 
have shown that in the Milky Way and external galaxies observed and real number of 
stationary sources may differ by 2 to 3 orders of magnitude. Given this uncertainty, our 
estimate of -/Vsss (Table |2J) does not contradict observations. 



5 Discussion 



Occurrence rate of accumulations of Mch by white dwarfs in semi-detached systems 
in our models does not exceed 0.22 x 10 _3 yr _1 (Table. 12 model 22). Even if we assume 
that the systems which in our model experience accretion-induced collapses actually 
explode as SNIa, the total rate of events does not exceed ~ 0.4 x 10 -3 yr" 1 (model 26); 
the latter model requires extremely high a ce and a kind of "optimization" on ath ! This 
model occurrence rate comprises only 20% of the lower limit of the observational estimate 
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of the occurrence rate of SNIa in our Galaxy. The model estimate is even somewhat 
overestimated (by 10 - 20%) since, even if the mass transfer rate is within the limits that 
correspond to the stationary burning of the hydrogen, some of the accreted matter has 
to be lost via stellar wind from a hot white dwarf. Nevertheless, the fraction of potential 
SN la that can occur in semi-detached systems is not negligible. A critical observational 
evidence for explosions of SNIa in low-mass semi-detached systems would be provided 
by discovery of hydrogen in the spectra of SN la. The main s ource of the latter sho uld be 



the matter stripped from the companion. Computations of iMarietta et alJ (|200Df ) have 



shown that for the systems under consideration the mass of stripped matter may be 
~ 0.15 Mq. Hydrogen may be observed both in very early and very late optical sp ectra, 
in radio, and in X-rays ( Eck et al. . 19951 : Marietta et al. . 2000l : Lentz et al. . 20021 ) . An 
unambiguous evidence in the favor of single degenerate scenario would be provided by 
discovery of the companion of exploding dwarf. This companion, due to the interaction 
with the envelope of SN, may have anomalously high luminosit y for the first ~ 10 3 — 10 4 
yr after explosion ( Marietta et al. . 2000l : Podsiadlowskil. l2003h . Spatial velocity of the 
companion after explosion may be up to ~ 450 km/s (jCanal et all . 200 ll ). 

For the estimate of the real role of the channel of semi-detached systems a crucial 
role may be played by an evaluation of the real number of the supersoft X-ray sources 

in the Milky Way and external galaxies. 

Recently discovered SNIa 2002ic dHamuv et all 1200.4 with hydrogen lines in the 
spectra is, most probably associated with a massive AGB object and may represent either 
an explosion of a single star (s o called SN 1.5 ) or an explosion of a degene r ate co mponent 
in a symbiotic binary system (jrlamuv et all 1200.4 Chugai fc Yungelsonl . 200.4 . 

Our estimates of the occurrence rate of accumulations of Mck by degenerate dwarfs 
seve ral times lower t han t h e estimates by Hachisu et alJ ( 199fil ) ; iLi fc van den Heuvell 



arc 



( 19971 ): lHachisu et ail (|l999h : lHan fc Podsiadlowskil (1200.4 who have found that the 



model estimates may be more comparable to the observational estimates. Hachisu et al. 
and Li and van den Heuvel have also used Eq. © for the estimate of the formation rate 
of binary systems, but did not carry out a detailed population synthesis study. Because 
of nonlinear variation of masses of components and semi- major axes of orbits in the 
course of evolution, simple replacement of numerical integration of the function © over 
the space of initial parameters by multiplication of the finite intervals of -/Wio, logao, 
and qo results in a significant overestimation of the occurrence rate of potential SNIa. 
Another reason for the difference with the results of other studies may be still uncertain 
efficiency of the matter accumulation by accreting dwarfs. 

For the "standard" assumption of a ce =l the same population synthesis code gives 
2.1 x 10~ 3 yr _1 for the rate of mergers of close double degenerates with total mass 
exceeding Mch- For the rate of edge-lit detonations the rate is 4.6 x 10~ 3 yr _1 . These 
estimates are comparable to the observational estimate of the occurrence rate of SN la 
in our Galaxy. 

The age of semi-detached systems considered in the present p aper at the moment 
of SN Ia explosion does not exceed ~ 2.5 Gyr (see for discussion Yungelson Sz Liyiol . 
1998). The age of SN initiated by edge-lit helium detonations does not exceed ~ 1 Gyr 
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( Yungelson fc Livid* . I2OO0I : Regos et all Eo 03) 2 . Hence, semi-detached systems cannot 



produce SN la in old populations (elliptical galaxies) . The merger of double-degenerates 
is apparently the only mechanism that is able to produce (super) Chandrasekhar mass 
objects in the populations of any age. 

The existence of close pairs of white dwarfs with total mass close to the Chan- 
drasekhar _jme_and^_ojj3ital peno^ enough for merger in Hubble time is out of 
question dNaniwotzki et all 120021 . l2003alHl. The r e lative number of such pairs is close 



to the theoretically expected one (jNelemans et all 120 01): 3 out of ~ 1000 field white 



dwarfs with V < 16™ 5 studied for binarity. However, it is still not clear whether the 
merger results in the carbon ignition in the center of the dwarf and SNIa. 

In the merger process, the less massive object turns into a disk that surrounds the 
"core". Further evolution of the system is determined by the rate of accretion from 
the disk M a . The rate of accretion depends on the time scale of angular momentum 
transfer in the disk T a . The latter may be of the order of 10 9 yr, if it is determined 
by the viscosity of degenerate electrons only. However, due to the turbulence in the 
transition layer betwee n the disk and the co r e, T„. is, most probably, by several orders 



of magnitude shorter ([Mochkovitch &: Livid . |199(J) and is comparable to the limiting 



Eddington accretion rate (M Edd ~ 2 x 10 5 M yr 1 for a 1M & white dwarf). 

The results of a ccretion computat ions are ambiguous even in ID case. 
Nomoto fc Ibenl dl98fJ ): iKawai et all (|1987l ) have shown that the central ignition of car- 



bon and SNIa after accumulation of Men is possible for M a (0.1 — 0.2) M Edd only. For 
larger M a , carbon ignites in the layer of accumulated matter and propagates inward. 
Ignition of carbon initiates the "competition" of two processes: propagation of carbon 
burning and matter accretion from the disk. This allows to get some simple estimates. 
The maximum rate of the carbon burning, (with account for the fraction a v of released 
energy that is taken away by neutrinos) is limited by the Eddington luminosity. Accre- 
tion rate also does not exceed the limit set by LEdd- The dwarf is able to accumulate 
Mch before the burning propagates to the center, if 

M C h ~ M a0 e g 

M Ch (1 - a u )X c e c ' 1 ' 

Here M a Q is the initial mass of accretor, X c ~ 0.5 is carbon abundance in the outer 
layers typical for massive dwarfs, ec ~ 5 x 10 17 ergg -1 is the energy released by burning 
of lg of carbon, and e g ~ 5 x 10 16 ergg -1 is the energy released by a ccretion of 1 g of 
the m atter from the disk. Since in the carbon burning layer a v ~ 0.4 ( Saio fc Nomotol . 



1998h , the estimate (J7J) shows that a dwarf with M a o = 1 M is able to accumulate M, 



Ch 

before carbon burns completely. Thus, it is possible that ignition of carbon in the outer 
layers of the dwarf does not prevent central explosion and SNIa. Only for a u ^ 0.9 
complete exhaustion of carbon results in formation of an oxygen-neon-magnesium white 
dwarf that, after accumulation of Mq^, "quietly", without SN, collapses into neutron 
star. 

2 It is possible th at the lifting effect due to rotation of the dwarf may prevent helium detonation 
ganger et all 1200.^ . 
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On the other hand, IPunina-Barkovskava fc Imsheniiikl |2Q03) have shown that, for 
the accretion rate of 5 x 10~ 7 M^yr" 1 and initial mass of CO white dwarfs (1.33 - 
1.34) Mq, carbon burning that starts at the edge of the dwarf and occurs in the regime 
of delayed detonation is possible. A SN la may result in this case. 

The situation may change if rotation is taken into account. iPiersanti et al. (2003) 



have shown that, for certain combinations of accretion rate and spin velocity of the 
white dwarf, rapidly rotating white dwarf deforms by adopting an elliptical shape. The 
anisotropic mass distribution induces the loss of rotational energy and angular mo- 
mentum via gravitational wave radiation. The white dwarf contracts and achieves the 
conditions suitable for explosive carbon burning at the center. However, adequate 3D 
computations are necessary for definite conclusions. 

6 Conclusion 

We have carried out a detailed study of the occurrence rate of accumulation of 
Mch by CO- and ONe-white dwarf components of semi-detached binaries. We have 
considered the dependence of model results on the parameter of common envelopes 
and conditions for formation of the common envelopes in the intervening stages of the 
evolution depending on assumptions on the mixing of He-rich accreted matter. We have 
shown that even for the most favorable combination of the parameters accumulation 
of Men by white dwarfs in the Galactic disk occurs at the rate that does not exceed 
~ 2 x 10 yr . Thus, in our model, this channel for formation of progenitors of 
potential SNIa is not able to produce more than ~ 10% of all SNIa in our Galaxy. An 
alternative channel for formation of presupernovae - the merger of double-degenerates 
is able to provide the necessary occurrence rate, but the details of the merger process 
itself and subsequent explosion mechanism need additional study. Among these details 
are, in particular, the structure and evolution of the "core+disk" configuration that 
is formed immediately after the merger, mass and angular momentum transfer in this 
configuration taking into account viscous processes and magnetic field, the character of 
burning that ensues after the merger, the role of neutrino losses and effects of rotation. 

In the semi-detached systems under consideration, "quite" accretion-induced col- 
lapses of ONe-white dwarfs may occur with the rate up to ~ 10~ 4 yr -1 . This rate, in 
principle, does not contradict the limits on the rate of AIC's imposed by ^^considera- 
tions of the production of r-process elements in the Galaxy dl-rvei- el all Il99flh . 



The authors are grateful to P. Kahabka for discussion of the properties of supersoft 
X-ray sources. This study was partially supported by RFBR grant no. 03-02-16254 and 
Federal Science and Technology Program "Astronomy". 

Appendix 

Let us make order of magnitude estimates of the time scales of several processes 
involved in t he accretion of He and compare them with the time scale of mass exchange 
t ex ~ 10 5 yr dlben fc Titukovl . Il98.4 iBitzaraki etaP . l2002h due to expansion of helium 



remnant of the primary component of initial system. It is natural to assume that, due 
to shear instability, the layer of accreted matter is turbulent. Let us assume that the 
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typical size of helium cells is equal to the local pressure scale height H p in the envelope 
of the star. 

The time scale of the radiative cooling of a turbulent bubble is 

Hp Hp Kp 

CI C 

where c is light velocity, I - mean free path of quanta, k - opacity, p - density. Or, in 
solar units, 

(R \ 4 / M \ 2 
m.) r2yr ' (A - 2) 

where R and M are radius and mass of the star, T - temperature. For a 2 M@ star 
tcool ~ 10 3 yr. It is evident that t coo i <C t ex . Hence, the temperature of a sinking 
element with the size equal to the local pressure scale height is practically equal to the 
temperature of the ambient matter. 

For the equal pressure, the equal temperatures result in the difference of the density 
of the bubble and the ambient matter and generation of the Archimedes force. Equating 
the force of dynamical friction irH^pv 2 and the Archimedes force (4ir/3)HpgAp, one may 
find the velocity of sinking: v ~ ct\/ Ap/ p, where ct is the local sound speed. Then the 
time of crossing the radius of the star by a turbulent element is tjf rs (R/ct) ^/ApJ/j,. 
Since for 2M Q star £// — 2hr, but the mass exchange time t ex ~ 10 5 yr, one obtains 
that the mixing occurs practically instantaneously, if relative difference of the molecular 
weight of the convective cell and stellar envelope is 



^ >10- u . (A.3) 

A* V tex J 

Equation (A.3) is valid always. 

Further, let us estimate the rate of diffusion of helium out of convective cells. If 
in the process of sinking the cell does not change its identity, its typical radius varies 
as r = (po/ pY'^HpQ, where H p q w 10 8 cm is the pressure scale height and po i=a 1CT 8 
gem" 3 is the density in the photosphere. The time scale of thermal diffusion may be 
estimated as 

r2 u l~3kT .. A . 

r d Ki — — , where v a = W . (A.4) 

v a L a y p a 

Here p a - mass of He nuclei, v a - their velocity, l a = l/an a - mean free path of the 
nuclei, n = p/p a - concentration, a = I6e 4 / (kT ) 2 cross-section of the interaction of 
nuclei, e - the charge of an electron ( Land . ^974^ . Inserting po and H p q given above, 
one gets 

r d «10 19 |^ yr. (A.5) 

The time scale Td — 10 2 — 10 5 yr is also much shorter than the mass exchange scale 
for whole star, with a possible exception of the most outer layers of the envelope with 
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relative radius of (AR/R) « 0.01. If the turbulent fragmentation of He-enriched sinking 
cells occurs, has to be even shorter. 

In the detailed modeling of the disk accretion of helium one has also have in mind 
the speed-up of stellar axial rotation and circulation of the matter caused by rotation 
that facilitate the mixing as well. 

The circumstances considered above, lead us to conclusion that the "instantaneous" 
mixing of whole star is quite probable and deserves further study. 
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Table 1: Parameters of the evolutionary tracks for helium-accreting star with initial 
mass of 2.Mq. 

dM/dt - accretion rate, R m ax/Ro & n d Rmin/Ro ~ the ratio of maximum and 
minimum values of the radius of accreting star and initial radius Rq, respectively. 
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Table 2: Occurrence rate of potential SNIa and of accretion-induced collapses of white 
dwarfs in semi-detached binaries as a function of parameters of computations. 
Observational estimate of the occurrence rate of SNIa in the Galaxy is (4 ± 2) x 10~ 3 yr 
l|(Jappellaro & Turattot I20TTH) . For every combination of a ce and a t h three upper lines give 
occurrence rate of possible SNIa and AIC's in scenario I under different assumptions on the 
evolution of the system at the helium accretion stage: conservative mass-exchange, formation of 
common envelope, and threefold expansion of accretor. The fourth line gives the rate of events 
in systems evolving through scenario II. -/Vggg is the number of supersoft X-ray sources in given 
model (for scenarios I and II together). 
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Figure 1: Schematic representation of a scenario of evolution resulting in the formation 
of a potential SNIa in the case B of mass-exchange ("a helium-rich supersoft X-ray 
source channel", scenario I). 
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Figure 2: Schematic representation of a scenario of evolution resulting in the formation 
of a potential SN la in the case C of mass-exchange, scenario II. 
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Figure 3: Variation of parameters of accreting stars depending on accretion rate and 
mode of mixing. Upper panel - tracks of the stars in the HR diagram. Lower panel - 
dependence of radii of accreting stars on mass. Notation corresponds to the numbers of 
tracks in Table ^ The position of initial model is marked by an asterisk. 
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Figure 4: Upper panel: variation of the profile of helium concentration in the models 
of track 2 with time (including the stages after termination of accretion). The numbers 
indicate time that elapsed since beginning of accretion in 10 3 yr. Lower panel: variation 
of stellar radii in the models of tracks 2 and 5 (solid curve corresponds to the stage 
of accretion while the dotted curve - to the evolutionary stage after termination of 
accretion. The letters indicate the models with helium profiles shown in the upper 
panel. 21 




Figure 5: Upper panel: progenitors of systems in which white dwarf can accumulate Mch 
in the "initial mass of the primary - initial separation of components" diagram. The 
border between systems experiencing cases B and C of mass exchange (for M20 = 2 M@) 
is shown by dashed line. The maximum of gray scale corresponds to the systems that 
provide SNIa occurrence rate of 0.94 x 10~ 5 yr _1 . Lower panel: the relation between 
mass of the donors and accretors in the systems in which white dwarf attains Mch at 
the moment when the donor overfills its Roche lobe (beginning of the stages 7 and 5 in 
scenarios I and II, respectively). The maximum of gray scale corresponds to the systems 

that provide SNIa occurrence rate of 0.15 x 10 -4 yr -1 . The figure corresponds to model 
13 in Table El 
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